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Surface probes such as scanning tunneling microscopy (STM) have detected complex patterns at 
the nanoscale, indicative of electronic inhomogeneity, in a variety of high temperature superconduc- 
tors. In cuprates, the pattern formation is associated with the pseudogap phase, a precursor to the 
high temperature superconducting state. Symmetry breaking (i.e. from C4 to C2) in the form of 
electronic nematicity has recently been implicated as a unifying theme of the pseudogap phase,— how- 
ever the fundamental physics governing the nanoscale pattern formation has not yet been identified. 
Here we use universal cluster properties extracted from STM studies of cuprate superconductors in 
order to identify the fundamental physics controlling the complex pattern formation. We find that 
the pattern formation is set by a delicate balance between disorder and interactions, leading to a 
fractal nature of the cluster pattern. The method introduced here may be extended to a variety of 
surface probes, enabling the direct measurement of the dimension of the phenomenon being studied, 
in order to determine whether the phenomenon arises from the bulk of the material, or whether it 
is confined to the surface. 



While lanthanum-based cuprate superconductors dis- 
play striking evidence of rotational symmetry breaking 
in the pseudogap regime, sometimes even leading to 
translational symmetry breaking^, such issues have been 
harder to settle in the higher transition temperature com- 
pounds such as YBCO and BSCCO. Recent experimental 
progress on YBCO has shed new light, including nematic 
behavior in the pseudogap regime detected via Nernst 
effecli^, transport^, and neutron scattering^, as well as 
evidence of time-reversal symmetry breaking detected via 
neutron scattering^ and the Kerr effecij^. The detection 
of a glass of unidirectional domains via STM in NCCOC 
and Dy-Bi2212 has now been followed by the dramatic 
demonstration of rather large electron nematic domains 
in BSCCOi^ 

Because of the fragile nature of rotational symmetry 
breaking, certain classes of disorder forbid long-range ori- 
entational order in 2D systems, and significantly suppress 
the order in strongly layered systems.— It is therefore 
critical to understand the effects of both quenched disor- 
der and interactions between nematic domains in these 
systems. As shown in Ref. [^, the orientational degree of 
freedom of the electron nematic in the cuprates in the 
presence of quenched disorder maps to a disordered Ising 
model in the following way. When an electron nematic 
forms, there is a preferred orientation to the electronic 
degrees of freedom. We consider Cu-0 planes which are 
locally C4 symmetric, with an incipient clclctron nematic 
which breaks the local rotational symmetry of the host 
crystal from C4 to C2, leading to two possible nematic 
orientations. We coarse grain the system, and define a 
local nematic order parameter by tr = ±1, correspond- 
ing to the two allowed orientations. The tendency for 
neighboring nematic regions to align is modeled as a fer- 
romagnetic nearest-neighbor interaction. 

Material disorder in the form of, e.g., dopant atoms 
competes with the ferromagnetic coupling between lo- 
cal nematic directors. Upon coarse graining, there are 



two broad classes of disorder which present themselves at 
the order parameter level: local energy density disorder 
(which includes random Tc disorder), and random field 
disorder. Local density disorder may arise in the form of, 
e.g., random bond disorder, in which the strength of the 
ferromagnetic coupling varies from coarse-grained site to 
coarse-grained site in the system. In addition, the lo- 
cal amplitude of the nematic order parameter can vary 
spatiallji. In an order parameter description, this type of 
disorder may be subsumed into randomness in the bond 
strengths. The other class of disorder, random field dis- 
order, arises when the pattern of doping atoms breaks 
rotational symmetry, thus favoring one or the other ori- 
entation of the nematic director in that region. This type 
of disorder couples linearly to the nematic director. The 
purpose of this paper is to gain an understanding, at the 
order parameter level, of both the relative importance of 
interactions and disorder, and also to discern the type 
of disorder, thus arriving at an order parameter level de- 
scription of the fundamental physics controlling the inho- 
mogencous pattern formation observed at the nanoscale 
in these materials. We consider a general model encom- 
passing both classes of disorder: 
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Here, sets the overall strength of the in-plane fer- 
romagnetic coupling between nearest neighbor Ising ne- 
matic variables, is an in-plane ferromagnetic coupling be- 
tween nearest neighbor Ising nematic variables, and 
represents bond disorder in the coupling strength. 
represents the overall coupling strength between Ising 
variables in neighboring planes, and SJ^^ is bond disorder 
in the interplane coupling strength. The random field hi 
is chosen from a gaussian probability distribution ccn- 
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tcrcd about zero, with width A, which we cah the "ran- 
dom field strength". We have chosen to discretize the 
system on a cubic lattice for simplicity, a choice which 
does not affect the universal properties of the model as 
long as the mean value of the random bond couplings 
remains positive. 

The field h = /lint + /ioxt represents an orienting field 
which breaks rotational symmetry. The external contri- 
bution hext may be achieved by the application of, e.g., 
magnetic fields, uniaxial pressure, or high currents among 
others. In the systems we have analyzed, data was taken 
in the absence of applied fields. Another source of finite h 
may be internal crystal effects /lint, such as the chains in 
YBCO may present. Such issues do not arise in NCCOC 
or BSCCO, and it is appropriate in both NCCOC and 
BSCCO to set the internal field /lint = in our model. 

The order parameter m ~ '^i describes the de- 

gree of oricntational order in the system. In the limit 
of zero random field strength and no bond disorder, the 
model has a finite-temperature continuous phase tran- 
sition for dimension d > 2, from a disordered phase to 
an ordered electron nematic, i.e. with long range ori- 
entational order < m >^ 0. However, for any finite 
random field strength, long range orientational order is 
forbidden in two dimensions. In three dimensions, there 
is a finite critical disorder strength of the random field, 
Af^ = 2.27J||iii In strongly layered systems, the criti- 
cal disorder strength is finite, but significantly reduced 
from the 3D valueJ^ In contrast with random field mod- 
els, weak random bond disorder does not forbid nematic 
order in two dimensions; rather, the phase transition is 
governed by clean Ising model exponents. In layered and 
three dimensional systems, the phase transition of the 
random bond model is controlled by the disordered fixed 
point "R". In the presence of both random bond and 
random field disorder, the universality class is that of 
the random field model. 



I. UNIVERSAL SPATIAL PROPERTIES OF 
CLUSTERS 

Fig. [T] shows our extraction of locally oriented (ne- 
matic) domains from scanning tunneling microscopy data 
on Bi2Sr2Dyo.2Cao.8Cu208+i (Dy-Bi2212) from Fig. S3 
of Refi^, reported as an "R-map," where at each posi- 
tion f, R{r,V) = I{r,+V)/I{f,-V) is the ratio of the 
tunneling current at positive voltage to that at negative 
voltage<^"— Based on the "R"-map, the presence of lo- 
cal, unidirectional domains of width 4ao were noted^, 
corresponding to the distance between "legs" of the 4ao- 
wide ladders, where Oq is the Cu-Cu distance within the 
Cu-0 planes. There is also a coexisting local density wave 
at Co, corresponding to the distance between "rungs" on 
the ladders. To effect the mapping to Ising nematic do- 
mains, we take a local spatial Fourier transform (FT) of 
the R-map in any given region, and then focus on the 
Oo local density wave. There are two "flavors" to this 




FIG. 1. Extraction of Ising domains from tunneling asymme- 
try maps in Dy-Bi2212, with Tc ~ 45K, taken at T = 4.2A'. 
Notice the large, percolating cluster (orange), with isolated 
flipped domains inside. 



density wave, either oriented along the a axis, or along 
the b axis, depending on the local orientation of the lad- 
ders. The relative weight of the two density waves can 
be discerned from comparing the Uo periodicity peaks in 
the local FT, comparing the weight of the Oo peak in the 
a direction to that in the b direction. We assign a lo- 
cal Ising variable based on this relative weight. We have 
checked that our results are insensitive to details such as 
the size of the FT window and the Ising lattice spacing<21 
Note that our assignment of the Ising cluster map is in- 
dependent of the local strength of the nematicity, which 
can be subsumed into randomness in the Ising couplings 
Jij, contributing to random bond disorder. 

The extracted cluster maps resemble cluster patterns 
observed in numerical studies of Eqn. [2p^ consistent with 
the idea of mapping disordered ncmatics to a disordered 
Ising modeli^ We show in Tables |T] and [TTl cluster prop- 
erties which can be extracted from the spatial configura- 
tions of clusters, as defined in the following references: i] 
is the spin-spin correlation function exponentii, r is the 
cluster size distribution exponenl^, and dc is the fractal 
dimension of the spanning cluster—. The theoretical val- 
ues available from the literature are shown for clean Ising 
models, random bond Ising models, and random field 
Ising models in both two and three dimensions ji^ii^i^ 
or derived from a combination of theoretical exponents 
available in the literature and scaling relations. (See the 
Appendix.) Notice how distinct the numbers arc from 
one universality class to the next. This implies that dif- 
ferent universality classes are sufficiently different that 
one can in principle use these methods to identify the 
universality class controlling the inhomogeneous pattern 
formation. In particular, because critical exponents are 
so sensitive to dimension, this method can in principle 
be used to allow any surface probe to directly measure 
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Clean Ising Model 


Random Bond Disorder 


Random Field Disorder 


T 


2.067-^'''' 


2.067--J^f^ 




dc 


1.125^i=i 


1.125iiiiii 




d-2 + r] 


0.25^ 


0.25ii 


1^ 



TABLE I. Theoretical Critical Cluster Exponents of Two Dimensional Ising Models. In two dimensions, random 
bond disorder is irrelevant, and the order to disorder phase transition is controlled by the clean Ising model critical point. For 
the random field Ising model in two dimensions, long range order is forbidden at finite disorder strength. However, there is an 
unstable fixed point at T=0 and zero disorder. The table reports values from the literature for this point, since it may aflect 
the scaling in some regimes. 





Clean Ising Model 


Random Bond Disorder 


Random Field Disorder 


r 


2.21— 


2.221— 


2.02±0.03i^'^ 


dc 


2.11-ti'iii 


1.98ii^ 


2.78 ± 0.05^ — 


d-l + j] 


1.04^ 


1.06i^ 


1.5±0.05^ 



TABLE II. Theoretical Critical Cluster Exponents of Three Dimensional Ising Models. 





Dy-Bi2212 at T=4.2K 


T 


1.71 ± 0.07 


dc 


1.56 ± 0.02 


d-2 + ri 


0.26 ± 0.03 



TABLE III. Critical Cluster Exponents from STM on 
the cuprate superconductor Dy-Bi2212. 
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methodsiii^"— from the statistical mechanics of 
disordered systems to identify the fundamental physics 
controlling the inliomogeneous pattern formation. Our 
results using published STM data,2i on Dy-Bi2212, a 
material which displays evidence of a disordered electron 
nematic at the surfacejii^ are shown in Table IIIIl Fig. [2] 
charts a direct comparison between our values extracted 
from data, and the theoretical fixed point exponents of 
disordered Ising models. Because of the strongly layered 
nature of Dy-Bi22I2, we expect there to be a regime of 
scaling which follows the 2D universality class, which 
then crosses over to the 3D universality class at longer 
length scales 



Cluster Size Distribution 



FIG. 2. Chart comparison of our results for Dy-Bi2212 com- 
pared to the theoretical models. Solid circles represent the 
values reported in the literature, summarized in Tables [I] and 
im Red circles represent r, blue circles represent dc, and green 
circles represent d — 2-\-r]. In comparing our results to the 3D 
models, we have assumed that at the surface, one observes a 
two-dimensional cross section of a 3D cluster, and we subtract 
1 from the fractal dimension dc. These values are represented 
by the open blue circles. Thick lines connecting circles repre- 
sent putative crossover of exponents from 2D to 3D behavior 
in a layered material. The horizontal lines in red, blue, and 
green are our results for the spatial characteristics of clusters 
in Dy-Bi2212 as reported in Fig. S3 of Ref. [H 

the dimension of the phenomenon studied, to determine 
whether the phenomenon is happening only on the sur- 
face of the material {d = 2), or whether it arises from the 
bulk of the material (d = 3). 

We apply quantitative cluster analysis 



The exponent r characterizes the distribution of clus- 
ter sizes. Fig. [T] shows the clusters identified from Dy- 
Bi2212. Here, a single cluster is a nearest-neighbor set of 
like-oriented Ising variables on the square lattice we have 
chosen to use. Fig. [3] shows the cluster size distribution, 
a histogram of cluster sizes. The results are noisier at 
large cluster size, due to the finite field of view (FOV). 
Fig m shows the same distribution with logarithmic bin- 
ning. Near a critical point, the cluster size distribution 
is of the form D{A) oc A~'^ . Using a straightforward 
fit to this form, we find r = I.7I ± 0.07. There is one 
large spanning cluster, represented by the last point in 
Fig. 131 Near criticality, scaling can be observed over a 
finite range. Note that in this case, the spanning clus- 
ter is also in the regime of scaling, and almost 4 decades 
of scaling are present in this cluster property, and this 
extraction of the exponent r can be considered quite re- 
liable. 
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FIG. 3. Raw cluster size distribution. 
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FIG. 4. Gluster Size Distribution after logarithmic binning, 
used to calculate the critical exponent r. 



B. Fractal Dimension of Clusters 

The fractal dimension of clusters relates the surface 
area of each cluster to its volume. Because STM is a sur- 
face probe, the available information is two-dimensional, 
thus the fractal dimension of clusters relates their perime- 
ter p to the area a of each cluster: p oc af-'^l'^w^ Using 
d = 2 because only a two dimensional cross section of 
the cluster properties is available, a straightforward fit 
gives dc = 1.56 ± 0.02. It is evident from Fig. [5] that 
the spanning cluster is also in the scaling regime for this 
measure, leading to almost 4 decades of scaling for the 
fractal dimension, and as with r, this extraction of the 
exponent dc can also be considered quite reliable. 



C. Spatial Correlations 

Within the Ising description, the spin-spin correlation 
function of the Ising pseudospin variables is described by 
Gir) (X \r\-(d.-'^+^) . Within the field of view available, 
we have averaged over all sites to obtain the spin-spin 
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FIG. 5. The fractal dimension of clusters, dc, relates the 
perimeter of clusters to their area. (See text for details.) 
Logarithmic binning has been used. 
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FIG. 6. Logarithmically binned data for calculating the crit- 
ical exponent rj. 

correlation function plotted in Fig. [51 At short distances, 
a weak region of scaling is evident, which extends roughly 
half of a decade. With only a half decade of scaling, this 
measure should be considered unreliable, especially since 
the presence of an as-yet undetermined scaling function 
may skew results at short distances. Nevertheless, for 
completeness, we note that the value of d — 2 -I- 77 obtained 
from this simple fitting procedure is d— 2+77 = 0.26±0.03. 



II. DISCUSSION 

The value of the cluster size distribution exponent r 
displayed by the clusters in Dy-Bi2212 lies between the 
2D and 3D values of the RFIM, indicating that the cluster 
size distribution is consistent with the scaling of a layered 
RFIM. By contrast, the value of t is 17% away from the 
clean and random bond values denoted by thick red lines 
in Fig. [2] Furthermore, the data show almost 4 decades of 
scaling for r indicating the value of t extracted is highly 
reliable. In comparing the fractal dimension of clusters dc 
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to 3D models, we make the ansatz that wc have analyzed 
a 2D cross section of a 3D cluster, and subtract 1 from 
the expected value of dc- We see then that the value of 
dc extracted from experiment is close to the theoretical 
RFIM values (within 12%), whereas it is 39% away from 
the random bond and clean models. There are almost 4 
decades of scaling for the fractal dimension dc, making 
this a reliable extraction of this universal cluster prop- 
erty. The value of c? — 2 + 77 extracted from a simple fit 
to the spin-spin correlation function is inconsistent with 
the random field universality classes, but it is consistent 
with Ising models, either clean or random bond mod- 
els close to two dimensions. However, this is the least 
reliable exponent extracted, as it displays less than one 
decade of scaling. The most reliable properties we have 
extracted are r and dc with almost 4 decades of scaling 
including the spanning cluster. (We note that even with- 
out the spanning cluster, over 2 decades of scaling are 
present in both of these quantities.) Overall, there is a 
striking similarity between the spatial properties of clus- 
ters observed in STM and universal cluster properties of 
disordered Ising models in 2D and 3D, particularly of the 
random field type. 

The presence of power law scaling in these quantities 
suggests that the pattern formation in Dy-Bi2212 is near 
criticality. Because the critical region is large in the 
RFIM it is possible to observe the scaling in a broad 
range of parameters near the critical point 1^ For exam- 
ple, two decades of scaling appear in the avalanche size 
distribution of the nonequilibrium 3D RFIM even at a 
disorder strength of A = 4 J, which is 85% larger than the 
nonequilibrium critical disorder strength = 2.16Ji^ 
The presence of power law scaling is also indicative of the 
fractal nature of this pattern formation, suggesting a po- 
tential connection to recent work on the fractal nature of 
doping in some high temperature superconductorspSS, 

In a true 2D condensed matter system {e.g. , with a fi- 
nite disorder concentration), random field disorder would 
forbid long-range order of an electron nematic. However, 
in any finite size 2D system, the system can appear to 
be ordered whenever the correlation length ^ is larger 
than the system size observed. Although the material we 
have considered is strongly layered, the layered RFIM ul- 
timately flows to the 3D universality class which has 
a finite critical disorder strength. Because of the weak 
coupling between planes, the critical disorder strength is 
suppressed from the 3D limit, i.e. Rc{J± < J\\ < Rc^)- 
Based on this, coupled with the fact that long-range- 
ordered nematicity has not been experimentally detected 
in BSCCO via a bulk probe, it is reasonable to expect 
the disorder strength in this material to be beyond the 
critical disorder strength. However, there is a spanning 
cluster in Fig. [U so the system cannot be too far above 
the critical disorder strength, since deep in the disordered 
phase, clusters are small. A spanning cluster was also ob- 
served in Ref. [H. We thus conclude that this material is 
in the intermediate disorder regimer^ where the disor- 
der strength is small within a plane, but strong between 



planes, J|| ^ A >> J±. Although long-range nematic 
order is not possible in this regime, clusters within each 
plane can have a long correlation length. 



III. CONCLUSIONS 

In conclusion, we have used quantitative cluster meth- 
ods from the statistical mechanics of disordered systems 
to identify the fundamental physics controlling the com- 
plex pattern formation revealed by STM in Dy-Bi2212. 
We find that the patterns are scale-invariant and fractal, 
indicating that Dy-Bi2212 is near a critical point, due 
to the interplay of local rotational symmetry breaking 
(electronic nematicity) with quenched disorder. We fur- 
thermore identify this critical point as being in the Ising 
universality class, with random field disorder the domi- 
nant type of disorder. Wc find that the disorder strength 
is in the intermediate regime, with disorder being a weak 
effect within each plane (leading to large clusters within 
a plane), while disorder is a strong effect in the interplane 
direction (leading to a lack of true long range order), in- 
dicating the pattern formation may best be described by 
a layered random field Ising model. Furthermore, this 
method may be extended to other materials and surface 
probes. For spatial patterns that may be tuned through 
or sufficiently near a critical point, this allows the de- 
termination of the dimension of the phenomenon being 
studied, enabling surface probes to distinguish whether 
the pattern formation is merely on the surface, or extends 
throughout the bulk of the material. 
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APPENDIX 

A. Theoretical Values of Critical Exponents 

The exact critical exponents we require in order to 
compare to known theoretical results of Ising models are 
not necessarily directly in the literature. However, scal- 
ing relations allow us to infer values for the cluster size 
distribution exponent, t, and the fractal dimension of 
clusters, dc. To infer dc, we employ the scaling relationii: 
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To infer r we used the scaling relationi^ 

2-a ^ 



P + 1 



(3) 



along with the Rushbrooke identitj*i^ to eliminate 7 in 
favor of a: 

a + 2/3 + 7 = 2^7=:2-2/3-Q; (4) 
to get the final expression: 

2 - a 



2- f3-a 



1 



(5) 



because /? is so small this equation is relatively insensitive 
to a 

Some of these exponents were extracted directly from 
the literature, in order to avoid confusion an asterisk will 
be used to denote that a particular value was inferred via 
relations rather than directly reported in the literature. 



Exponent 


2D Clean Ising Model 


3D Clean Ising Model 


d-2 + ri 
a 

P 

V 

T 

dc 


1/4^ 

1/sM 
■^33 

31/15* 

9/8* 


1.04iii 

0.1118(6)i^ 
0.326(4>i^ 
0.6294(2)i^ 

2.21* 

2.11* 



TABLE IV. Theoretical values of critical exponents in 2D and 
3D clean Ising models. 



Exponent 


2D RFIM 


3D RFIM 


d - 2 + 


1— 


0.5 ± 1.03^ 


a 




0.63^ 


P 

V 
T 


1.6i^ 


0.017(5)ii 
1.37(9)11 

2.01* 


dc 


1.9i& 


2.74* 



TABLE VI. Theoretical values of critical exponents in 2D and 
3D random field Ising models. 



Harris criterion considerations 



Exponent 


3D RBIM "R" point 


d-2 + ri 


1.029i^ 


a 


-0.042i 


P 


0.35^ 


7 


1.34^ 


V 


0.68^ 


T 


2.21* 


dc 


1.99* 



TABLE V. Theoretical values of critical exponents in 3D Ising 
models with weak random bond disorder. Note that transition 
in the 2D random bond Ising model is controlled by the clean 
fixed point. 



The Harris criterion states that local energy density 
disorder (such as random Tc disorder and random bond 
disorder) is irrelevant if dv > 2, where the exponents 
refer to the clean model. In the presence of hyperscaling 
(obeyed by the clean and random bond cases), dv = 2 — a 
implies that the Harris criterion reduces to a < 0. In 
the 3D clean Ising model, a k, 0.1—, and randomness 
is relevant. In the 2D clean Ising model, a = and 
such randomness is marginal. For the 3D case with weak 
bond disorder, there is a disordered fixed point with new 
exponents. In the 2D case with weak bond disorder, it 
has been shown that the system flows toward the clean 
model-^i^. 
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